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Abstract 

A simple method to calculate the total annihilation cross section of massive cold relics 
is presented. Exact expressions for a re i, b Te i, and c re i of av re \ = a re \ + b ve \v^ el + c ie ivf el are 
obtained from simple derivatives of the squared annihilation amplitude. Three examples are 
given. In one example, the computed annihilation cross section of massive Dirac neutrinos 
confirms recent correct calculations. 
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1 Introduction 



The number density of cold relics from the early universe depends sensitively on the anni- 
hilation cross section of such particles. The equation governing the evolution of the number 
density is the Boltzmann rate equation 



where n is the number density, H is the Hubble constant, no is the number density in thermal 
equilibrium and < av m \ > is the thermal averaged annihilation cross section. The thermal 
averaged cross section [] is approximately |2j 



where x = T/m (the cold relic mass is denoted as m). The freeze-out temperature is defined 
as the temperature at which the expansion rate of a co-moving volume becomes larger than 
the rate of annihilation. After freeze-out the relic particle decouples from the background, 
and the subsequent relic density is determined solely by the expansion rate of the universe. 

As has been shown by many authors the freeze-out temperature can be very small so as 
to justify using the Boltzmann distribution to calculate the thermal average. Most notably 
this is true in the case of stable neutralinos where Xf ~ 1/20 [[J. This low decoupling energy 
leads one to consider a low velocity expansion of the annihilation cross section: 



Obtaining the thermal averaged annihilation cross section then becomes merely a calculation 

of < >. However, when the relic particles annihilate into a narrow resonance special 

care must be taken [?]]. Although the low velocity expansion breaks down under these 

conditions p, |J, in most cases approximating the annihilation cross section with Eq. (|3]) is 

extremely useful. A crucial task remaining is the determination of a re \, b Ie \, and c Te \. 

In this letter I will demonstrate a simple technique to calculate the low velocity expansion 

of the annihilation cross section of 2 — > 2 processes with arbitrary final state masses. The 
2 For a precise treatment of the thermal average see refs. || and || 






(2) 



av Te \ = a re i + b re \v rel + c re iw rel H . 



(3) 
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technique substantially increases the facility of annihilation cross section calculations. Instead 
of performing integrals to obtain av ve \ and then expanding in v Te i, I expand in v re \ and then 
do the remaining integrals over simple trigonometric functions analytically to arrive at the 
final answer. This technique is made possible by the fact that the Mandelstam variables 
have no angular dependence in the zero velocity limit. Therefore, the differential cross 
section can have no angular dependence in the zero velocity limit, and can be expanded 
in a Taylor series with each power of velocity having coefficients of simple trigonometric 
functions. These simple trigonometric functions are then integrated out painlessly to obtain 
the total annihilation cross section. 



arbitrary final state masses is presented. The incoming momenta are labelled p\ and P2, 
both with mass m; the outgoing momenta are labelled p^ and p<± with masses m% and 
respectively; and, v cm is the center of mass velocity of the annihilating particles (v cm = v re \/2). 
The differential cross section is 



2 Expansion of the Annihilation Cross Section 



In this section the low velocity expansion of the 2 



2 annihilation cross section with 



da \M\ 2 |p 3 



(4) 



First, the following definitions are useful Q 




(5) 



Since 



Pil = 




mv, 



cm 



and 





the differential cross section can be written as 




(7) 



3 Recall that u = J2 m i~ s ~t ■ 
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The low velocity expansion of the Mandelstam variables is straightforward: 

s = s + s 2 v 2 cm + s 4 v* m + 0(v 6 cm ) 



where 



t = t + t 1 COS 0v cm + t 2 V 2 cra + t 3 COS Ovl m + t 4 V* m + O(^cm) 



So — s 2 = S4 = 4m 2 
—2m 2 + 777.3 + m l 

to = 2 



* 2 = -2m 2 



2m 2 (4m 2 — m§ — m 2 ) 



t 4 = -2m 2 



with 



/i 2 = ^/ (4m 2 — (m 3 — m 4 ) 2 )(4m 2 — (m 3 + m 4 ) 2 ). 
Next, expand K(s) in a Taylor series: 

K(s) =K + K 2 v 2 cm + K A vi m + 0{vtJ 



where 



K 
K 2 



2567rm 4 



2567rm 4 



' 3 »m 
2567rm 4 \8 + 74 



3 2m 2 ^ 2 \ 
2 + ^Tj 

1 m 2 fil 2m 4 /4" 



^1 



with /x| defined to be 



= 8m 2 — (m 3 — m 4 ) 2 — (m 3 + m 4 ) 2 . 

Expanding J(s,t)/4ir about s = s and t — t yields 
J(s,t) 



An 



^Ko + ^( S - S o) + ^(t-t ) 



+ 



i a 2 Jo 

2 9s 2 



(s - s ) 2 + 



l d 2 J 
2 <% 2 



9 2 J 



where Jo = J(s ,t ). The key point here is that Jo and d m+n Jo/ds m dt n do not depend 
on cos#; only (t — to) n can depend on cos 8. Therefore, to fourth order in velocity only the 
functions cos n 9 with n = . . . 4 need to be integrated. Of course / dfl cos 9 = J dQ cos 3 9 = 
which forces all terms in av rc \ which are odd powers of v re i to be zero. The result of integrating 
J(s, t)/4ir in the low velocity limit is 

dQ-IipH = J + jrf + j iV * + (v e c J (23) 



where 



47T 



Jo = J(s ,to), (24) 
d.J d.J 1 2 d 2 J 

J2 = 52 97 + ^-W + B^-dF' (25) 

dJ 1 2 d 2 J 1 d 2 J 

Ji — t^—, — h —t 2 „ „ + -t\t^- 



dt 2 z dt 2 3 1 ° dt 2 

1. ^ J , 1 +4 <9 4 J n 9 Jo 



Jp . 1 , 2 g 3 JQ . 1 2 ^Jo 

! 9 S ^ + 6 S2 1 9 S ^ 2 + 2 S2 9s 2 



+^2^ + ,s 2 t 2 ° + -s 2 , f (26) 



Expanding out all the terms in Eq. (0) the total annihilation cross section to fourth order 
in velocity is 

<?v cm = a cm + b cm v 2 m + c cra v* m (27) 

where 

Ocm = -«rcl = Jo-^0 (28) 

b cm = 2b ml = J K 2 -^J K + J 2 K (29) 

11 1 

c C m = 8c re i = J K 4 - ~ J o K 2 ~ sJqKq + J 2 K 2 - -J 2 K + J 4 K . (30) 

Thus, in principle the coefficients a cm , b cm , and c cm are easily calculated once the squared 
amplitude for the annihilation process is known. Although the derivatives of J(s,t) can be 
quite complicated, they are nevertheless always tractable. One is now faced only with the 
algebraic problem of taking derivatives. 
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3 Applications of the General Formula 



3.1 Annihilation of Massive Neutrinos 



Limits on the mass of neutrinos can be inferred from requiring that their relic density not 
overdose the universe. These bounds depend sensitively on their annihilation cross section. 
As a check to the above formalism, I calculate the annihilation cross section of both Majorana 
neutrinos and Dirac neutrinos into fermions and compare with previous calculations in the 
literature. In the case of Dirac neutrinos I confirm a recent result which disagrees with 
earlier literature. 

The interaction Lagrangian of Majorana neutrino annihilations into fermions is 



G F _ 

which corresponds to the following squared amplitude: 



(31) 



J{s,t) = AG 2 F {C\ + C^nf + 2m)-Am 2 s + s 2 - Am 2 t - Arnit + 2st + 2t 2 } (32) 



+ 16G : 



(5C 2 A - C 2 v )m 2 m) - C\m)i 



where m is the mass of the neutrino and m/ is the mass of the fermion. 

Using the technique developed in the previous section, the values of a Te \ and fe re i are easily 
calculated to give the following annihilation cross section to 0(v 2 el ): 



<JV re i 



Orel + fyrel^rel 

2G% 



-m 



71 



l-m 2 /mH{C 2 +Ciy-f 



C\ mi 
H — -— + 
2 m 2 



C 2 v - l ^C 2 A + 



6 

r 2 



m 2 12 



(33) 
and 



4(1 — m 2 /m 2 ) 

The above equation is identical to the annihilation cross section calculated by ref. 
others. 

The Dirac neutrino annihilation cross section is calculated similarly. The interaction 
Lagrangian is 



G 



F _ 



V2 



vY(l- l5 )vf^(C v -C Al5 )f 



(34) 



which gives the following squared amplitude: 

J(s,t) = 4G 2 F {C 2 A + C 2 ){m' l + m 4 f -2m 2 t~2m 2 f t + t 2 ) 



+2G 2 F (C V - C A ) 2 (s 2 + 2st - 2m 2 s) 



+8G 2 F (2C 2 A m 2 mj - C 2 A mjs + C A C v mjs). 



The low velocity expansion of av Te \ then becomes 
— m"\l 1 — m f 



av re i 



?mVl-m 2 f /m 2 (C 2 + C 2 A , 



24 ^ m 2 1 — m 2 /m 2 



m 



1 

2' 



1 + 



-rcl 



1 + 



1 



(35) 



(36) 



This result differs from refs. 
and Rothstein PI. 



1 — m 2 /m 2 

but confirms the recently published calculation of Dolgov 



3.2 Neutralino s— channel Annihilation 



Ellis et al. || calculated the thermal averaged cross section to 0(x) for the s-channel Z boson 
contribution to neutralinos annihilating into fermions. As a final example to demonstrate 
the technique developed in this paper, I calculate the annihilation cross section and compare. 
Following the conventions of ref. M the interaction Lagrangian is 



C = £(X7 M 75X) 
/ 



(g/g ~ %(jvl m2 z) m \ 

s — m 2 z + iTzmz 



[ff{A z P L +B z P R )f]. 



(37) 



From this Lagrangian the squared amplitude is calculated, and the values of J , J2, K , 
and K 2 are determined. Using Eq. (0) I find the coefficients a x and b x of < av re \ >= a x + b x x 
to be 

/ ... 9 A 

m . 



E 



1- 



??? 



/ 



ml 



27r V m 2 (m| — Am 2 ) 2 + r|m| 



m z 



(3? 



4 

m| 



2^ 



-Am 2 



m 2 (m| — Am 2 ) 2 + Y z m z 



Am 2 {A 2 z + B 2 z ) 



A 2 Z - 3A Z B Z + B z -3 



\m z J 



z — B z ) 2 



+ (~P f + 12 7/ + 5 f )m 2 {A z - B z ) 



( m V 
1-4 — 

\m z / 



(39) 



where 



rrif m 2 (m| — 4m 2 ) 3 

Pf = ™2 _ ^2 > 7/= C m | _ 4m 2)2 + r 2 m 2 ' 5 /=2' " 0) 



m" — m 



/ 



The value of a x obtained is identical to the expression for az found by Ellis et al. For b x I 
find that 

3 

bz = b x - 5 f a x = b x - -a x . (41) 

This discrepancy can be accounted for by different thermal averaging prescriptions. The 
thermal averaging procedure of ref. ||, which was used by ref. ||, requires an additional 
factor of 3a x /2 to be subtracted from b x . 

4 Conclusion 

Calculating annihilation cross sections accurately is extremely important since they, among 
other things, sensitively affect cosmological mass bounds of massive cold relics. Unfortu- 
nately, it is rare to find two independent calculations in the literature which agree with 
each other on a particular annihilation cross section. The technique outlined in this letter 
reduces substantially the possibilities of errors, and also greatly reduces the amount of effort 
required to calculate annihilation cross sections in the low velocity limit. Furthermore, using 
a symbolic manipulation program with this technique reduces even the most complicated 
processes to a few seconds of CPU time. 



Note Added: While completing this work I was made aware of a similar approach to this 
problem [|10|. 
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